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Differential measurements of charged particle azimuthal anisotropy are presented for lead-lead 
^ ! collisions at ^/£nn = 2.76 TeV with the ATLAS detector at the LHC, based on an integrated luminosity 
of approximately 8 /itr 1 . This anisotropy is characterized via a Fourier expansion of the distribution 
of charged particles in azimuthal angle relative to the reaction plane, with the coefficients v n denoting 
the magnitude of the anisotropy. Significant v 2 -v$ values are obtained as a function of transverse 
momentum (0.5 < p T < 20 GeV), pseudorapidity (|r/| < 2.5) and centrality using an event plane 
method. The v n values for n > 3 are found to vary weakly with both 77 and centrality and their p T 
dependencies are found to follow an approximate scaling relation, vn n (pr) oc v^ipr), except in the 
top 5% most central collisions. A Fourier analysis of the charged particle pair distribution in relative 
azimuthal angle (A0 = a - 4>b) is performed to extract the coefficients v n , n = (cosnA0). For pairs 
of charged particles with a large pseudorapidity gap (|Ar? = ?? a - •%{ > 2) and one particle with 
p T < 3 GeV, the ^2,2-^6,6 values are found to factorize as v nj „(p T ,p T ) « v n (p^)v n {p\) in central and 
mid-central events. Such factorization suggests that these values of ^2,2-^6,6 are primarily due to the 
response of the created matter to the fluctuations in the geometry of the initial state. A detailed study 
shows that the ui,i(p|,p T ) data are consistent with the combined contributions from a rapidity-even v\ 
and global momentum conservation. A two-component fit is used to extract the v\ contribution. The 
extracted v\ is observed to cross zero at p T ~ 10 GeV, reaches a maximum at 4-5 GeV with a value 
comparable to that for v 3 , and decreases at higher p T . 
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Differential measurements of charged particle azimuthal anisotropy are presented for lead-lead 
collisions at -^/snn = 2.76 TeV with the ATLAS detector at the LHC, based on an integrated 
luminosity of approximately 8 /_tb _1 . This anisotropy is characterized via a Fourier expansion of the 
distribution of charged particles in azimuthal angle relative to the reaction plane, with the coefficients 
v„ denoting the magnitude of the anisotropy. Significant «2-«6 values are obtained as a function of 
transverse momentum (0.5 < px < 20 GeV), pseudorapidity (\t]\ < 2.5) and centrality using an event 
plane method. The v n values for n > 3 are found to vary weakly with both 77 and centrality, and their 
Px dependencies are found to follow an approximate scaling relation, v„ (Pt) °t v 2 (px), except 
in the top 5% most central collisions. A Fourier analysis of the charged particle pair distribution in 
relative azimuthal angle (A(f> = 4> a — 4>h) is performed to extract the coefficients v n ,n — (cosnA</>). 
For pairs of charged particles with a large pseudorapidity gap (\Arj — r\ & — rfa\ > 2) and one 
particle withpT < 3 GeV, the «2,2-«6,6 values are found to factorize as «h,ti(Pt,Pt) ~ v n (p^)v n (p^) 
in central and mid-central events. Such factorization suggests that these values of t)2,2-t>6,6 are 
primarily due to the response of the created matter to the fluctuations in the geometry of the 
initial state. A detailed study shows that the t>i,i (Px,Px) data are consistent with the combined 
contributions from a rapidity-even vi and global momentum conservation. A two-component fit is 
used to extract the t)\ contribution. The extracted v\ is observed to cross zero at px ~ 1.0 GeV, 
reaches a maximum at 4-5 GeV with a value comparable to that for vz, and decreases at higher px- 
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I. INTRODUCTION 

The primary goal of high-energy heavy ion physics is 
to understand the properties of the hot and dense mat- 
ter created in nuclear collisions at facilities such as the 
Relativistic Heavy Ion Collider (RHIC) and the Large 
Hadron Collider (LHC). An important observable to- 
wards this goal is the azimuthal anisotropy of particle 
emission. At low px (;$ 3-4 GeV), this anisotropy re- 
sults from a pressure-driven anisotropic expansion of the 
created matter, with more particles emitted in the direc- 
tion of the largest pressure gradients [1]. At higher px, 
this anisotropy is understood to result from the path- 
length dependent energy loss of jets as they traverse the 
matter, with more particles emitted in the direction of 
smallest path- length [2]. These directions of maximum 
emission are strongly correlated, and the observed az- 
imuthal anisotropy is customarily expressed as a Fourier 
series in azimuthal angle <fi [3, 4]: 
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where px is the transverse momentum, 77 is the pseu- 
dorapidity, and v n and $„ represent the magnitude and 
direction of the n th -ordcr harmonic, respectively (see Sec- 
tion IV A). The n th -order harmonic has n-fold periodicity 
in azimuth, and the coefficients at low px are often given 
descriptive names, such as "directed flow" (v\), "elliptic 
flow" (v 2 ), or "triangular flow" (w 3 ). 
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In typical non-central heavy ion collisions where 
the nuclear overlap region has an "elliptic" shape 
(or quadrupolc asymmetry) on average, the azimuthal 
anisotropy is expected to be dominated by the «2 com- 
ponent [5-7]. However, it was recently pointed out that 
the positions of the nucleons in the overlap region can 
fluctuate to create matter distributions with additional 
shape components, such as dipole (n = 1) and sextupolc 
(n = 3) asymmetries [8-11]. Due to strong final-state in- 
teractions, manifested as either pressure or jet energy 
loss, these spatial asymmetries can be converted into 
final-state momentum anisotropics, leading to non-zero 
first-order and higher-order harmonic coefficients [11, 12]. 

The observation of large V2 for px < 3-4 GeV at the 
RHIC [13, 14] and LHC [15, 16] has led to the conclusion 
that the hot and dense medium behaves like a "perfect 
fluid" [14, 17, 18]. This is because the large V2 values re- 
quire hydrodynamic models [19-21] with a shear viscosity 
to entropy density ratio that is close to the conjectured 
lower bound of 1/4-7T [22, 23]. Precise determination of 
this ratio using only v-i data is limited by many model 
uncertainties [24]. Because the shear viscosity tends to 
dampen the harmonics, with more damping for larger 
n [11, 25, 26], measurements of harmonic coefficients be- 
yond V2 can provide stronger constraints for the shear 
viscosity of the medium. Extending these measurements 
to higher px is also valuable for discriminating between 
jet-quenching models, as high-px v n is sensitive to the 
path- length dependence of the jet energy loss [27-29]. 
These coefficients can also help to distinguish between 
different models of the initial geometry [30-33] , and pro- 
vide insights into the granularity of the initial state fluc- 
tuations [26, 34-36]. 



Another related observable for studying the properties 
of the medium is the correlation function between two 
particles in relative azimuthal angle A</> — <\> & — 0b and 
pseudorapidity A77 = r) a — 77b [37]. The distribution of 
pairs in Acj) can be expanded into a Fourier series: 
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where the coefficients v n ^ n are symmetric functions with 
respect to p T and p T . The harmonics defined in Eq. 1 
also contribute to this distribution: 
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where the global direction $„ drops out in the convolu- 
tion, and v„ is assumed to be independent of 77 (which is 
approximately true within \q\ < 2.5 at the LHC, see Sec- 
tion VA). Thus if the anisotropy is driven by collective 
expansion, v n ,n should factorize into the product of two 
single-particle harmonic coefficients [37]: 



Vn,n{PT,P%) = V n (jfe)v n (p\) 



(4) 



Such factorization may also be valid if the anisotropics 
of the two particles arc independently driven by collec- 
tive expansion and path-length dependent jet energy loss 
(both are associated with the same initial spatial asym- 
metries). This factorization relation has been used to 
calculate the single-particle v n [38-40]. On the other 
hand, autocorrelations induced by resonance decays or 
fragmentation of back-to-back jets, are expected to break 
the factorization. Therefore, Eq. 4 can be used to iden- 
tify the regions of p^ and p\ where correlations are dom- 
inated by effects controlled by the initial spatial asym- 
metries. 

The study of the structures of two-particle correlation 
in A77 and A</> has been the focus of major experimen- 
tal and theoretical efforts in the last decade. In typi- 
cal proton-proton collisions, where a medium is presum- 
ably not formed, the pair distributions arc dominated 
by strong correlation peaks at (A</3, A77) ~ (0, 0) and 
A(f) ~ 7r. These peaks reflect mainly autocorrelations 
among particles from fragmentation of back-to-back jets. 
In heavy ion collisions, additional structures have been 
observed for p^ < 3-4 GeV and large A77 at Ac/> ~ 
(known as the "ridge") [41, 42] and \A<f> — tt| — 1.1 
(known as the "double-hump") [37, 43]. These unex- 
pected structures have been interpreted as the response 
of the medium to the energy deposited by quenched 
jets [44, 45]. However, similar structures can also be 
generated by the flow harmonics, as they all contribute 
constructively at A<j) ~ but tend to cancel on the away- 
side according to Eq. 3 [8] . Therefore, a detailed compar- 
ison between the measured pair distribution (Eq. 2) and 
that expected from anisotropic flow (Eq. 3) can deter- 
mine whether the structures in two-particle correlations 
are a consequence of the so-called "jet-induced medium 



response" , or whether they are a consequence of a sum 
of the flow harmonics. 

The i>2 coefficient has been extensively studied at the 
RHIC [13, 28, 46-49] and LHC [15, 16]. Results for 
higher-order v n for n > 3 also became available re- 
cently [30, 31, 40]. In contrast, no experimental mea- 
surement of v 1 including systematic uncertainties exists 
at the LHC, although an estimate has recently been per- 
formed by a theoretical group [50] based on published 
ALICE data [39]. A primary complication for v\ mea- 
surements is global momentum conservation, which in- 
duces a significant dipole component [51, 52]. A "side- 
ward" deflection of colliding ions can also lead to a small 
rapidity-odd (i.e. changes sign crossing 77 = 0) dipole 
component [53, 54]. Therefore, the extraction of v\ val- 
ues associated with the initial dipole asymmetry requires 
careful separation of these contributions, which generally 
break the factorization relation given by Eq. 4. 

This paper presents comprehensive results for v\-vq 
over broad ranges of centrality, pseudorapidity and pt 
for charged particles in lead-lead (Pb-Pb) collisions at 
y/s^i = 2.76 TeV with the ATLAS detector at the LHC. 
The v n values are measured directly using an "event 
plane" (EP) method for n = 2-6, and are also derived 
from the v n>n measured using a two-particle correlation 
(2PC) method for n = 1-6. These detailed measure- 
ments provide new insights into the hydrodynamic pic- 
ture at low pt, the jet energy loss picture at high px, 
and the nature of the fluctuations in the initial geome- 
try. They also allow a detailed study of the factorization 
relation (Eq. 4) over broad ranges of centrality, A77, pfj-i 
and p\. Together, these measurements should shed light 
on the physics underlying the structures observed in two- 
particle correlation functions. 

The paper is organized as follows. Sections II and III 
give a brief overview of the ATLAS detector, trigger, and 
selection criteria for events and tracks. Section IV dis- 
cusses the details of the EP method and the 2PC method 
used to measure the v n . Section VA presents results for 
V2~vq from the EP method as a function of p^ , 77 and 
centrality. Section VB presents a detailed Fourier anal- 
ysis of the two-particle correlation functions to measure 
v n>n as a function of p T , p T , A77 and centrality, which 
are then used to calculate v 2—v q via the factorization re- 
lation (Eq. 4). These v n values are compared with those 
obtained from the EP method in Section VC, with a fo- 
cus on understanding the structures of the 2PC in terms 
of single-particle v n . Section V D presents results for V\ 
based on a two-component fit of the v\,\ data with a 
modified functional form of Eq. 4 that includes the con- 
tribution of global momentum conservation. Section VI 
gives a summary of the results and main observations. 



II. ATLAS DETECTOR AND TRIGGER 

The ATLAS detector [55] provides nearly full solid an- 
de coverage of the collision point with tracking detectors, 



calorimeters and muon chambers, well suited for mea- 
surements of azimuthal anisotropics over a large pseudo- 
rapidity range . This analysis primarily uses three sub- 
systems for v n measurement: the inner detector (ID), the 
barrel and endcap electromagnetic calorimeters (ECal) 
and the forward calorimeter (FCal). The ID is contained 
within the 2 T field of a superconducting solenoid mag- 
net, and measures the trajectories of charged particles 
in the pseudorapidity range |t;| < 2.5 and over the full 
azimuth. A charged particle passing through the ID typ- 
ically traverses three modules of the silicon pixel detec- 
tor (Pixel), four double-sided silicon strip modules of the 
semiconductor tracker (SCT) and, for \r]\ < 2, a transi- 
tion radiation tracker composed of straw tubes. The elec- 
tromagnetic energy measurement of the ECal is based on 
a liquid-argon sampling technology. The ECal covers the 
pseudorapidity range |rj| < 3.2, and is used as a reference 
detector in the event plane measurements. The FCal con- 
sists of three longitudinal sampling layers and extends the 
calorimeter coverage to \tj\ < 4.9. It uses tungsten and 
copper absorbers with liquid argon as the active medium, 
and has a total thickness of about 10 interaction lengths. 
The centrality measurement uses towers in all three lay- 
ers of the FCal, while the event plane measurements use 
towers in the first two layers of the FCal excluding those 
at the edge of the FCal 77 acceptance. These selection 
criteria are found to minimize the effect of fluctuations 
in the reaction plane measurement. 

The minimum-bias Level- 1 trigger used for this anal- 
ysis requires signals in two zero-degree calorimeters 
(ZDC), each positioned at 140 m from the collision point, 
detecting neutrons and photons with |?/| > 8.3, or ci- 
ther one of the two minimum-bias trigger scintillator 
(MBTS) counters, covering 2.1 < \r/\ < 3.9 on each side 
of the nominal interaction point. The ZDC Level-1 trig- 
ger thresholds on each side are set below the peak cor- 
responding to a single neutron, e.g. as produced from 
Coulomb dissociation of the lead ion [56]. A Lcvel-2 tim- 
ing requirement based on signals from each side of the 
MBTS is imposed to remove beam backgrounds. 



III. EVENT AND TRACK SELECTIONS 

This paper is based on approximately 8 ^b _1 of Pb-Pb 
data collected in 2010 at the LHC with a nucleon-nuclcon 
center-of-mass energy ^snn = 2.76 TcV. An offline event 
selection requires a reconstructed vertex and a time dif- 
ference | At | < 3 ns between the MBTS trigger counters 



1 ATLAS uses a right-handed coordinate system with its origin 
at the nominal interaction point (IP) in the center of the de- 
tector and the z-axis along the beam pipe. The x-axis points 
from the IP to the center of the LHC ring, and the y-axis points 
upward. Cylindrical coordinates (r, (p) are used in the trans- 
verse plane, (f> being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 9 as 
r) = -lntan(0/2). 



on either side of the interaction point to suppress non- 
collision backgrounds. A coincidence between the ZDCs 
at forward and backward pseudorapidity is required to 
reject a variety of background processes, while maintain- 
ing high efficiency for non-Coulomb processes. 

Events satisfying these conditions arc then required 
to have a reconstructed primary vertex within |z v t x | < 
150 mm of the nominal center of the ATLAS detector for 
the EP analysis. A more stringent vertex cut of 100 mm 
is required for the 2PC analysis, such that enough events 
can be found in the same z v tx bin for the event mixing 
procedure (see discussion in Section IV B). About 48 mil- 
lion and 43 million events pass the requirements for the 
EP and 2PC analysis, respectively. Pile-up probability is 
estimated to be at the 10~ 4 level and is therefore negli- 
gible. 

The Pb-Pb event centrality is characterized using the 
total transverse energy (J^-Et) deposited in the FCal 
over the pseudorapidity range 3.2<|?7|<4.9at the elec- 
tromagnetic energy scale. An analysis of this distribution 
after all trigger and event selections gives an estimate of 
the fraction of the sampled non- Coulomb inelastic cross- 
section to be 98 ± 2% [57]. This estimate is obtained 
from a shape analysis of the measured FCal J2 Et dis- 
tributions compared with a convolution of proton-proton 
data with a Monte Carlo Glauber calculation [58]. The 
FCal J2 Et distribution is then divided into a set of 5% 
or 10% percentile bins, together with a bin defined for 
the 1% most central events. The uncertainty associated 
with the centrality definition is evaluated by varying the 
effect of trigger and event selection inefficiencies as well 
as background rejection requirements in the most periph- 
eral FCal Y,Et interval [57]. 

Tracks are reconstructed within the full acceptance of 
the ID, requiring pr > 0.5 GcV and \r)\ < 2.5. To im- 
prove the reliability of track reconstruction in the high- 
multiplicity environment of heavy ion collisions, more 
stringent requirements on track quality, compared to 
those defined for proton-proton collisions [59], are used. 
At least nine hits in the silicon detectors (out of a typical 
value of 11) are required for each track, with no missing 
Pixel hits and not more than one missing SCT hit, in 
both cases where such hits are expected. In addition, 
the point of closest approach is required to be within 
1 mm of the primary vertex in both the transverse and 
longitudinal directions [16]. This selection is varied in the 
analysis to check the influence of both the acceptance and 
fake tracks. The tracking efficiency for charged particles 
is studied by comparing data to Monte Carlo calcula- 
tions based on the fflJING event generator [60] and a 
full GEANT4 [61] simulation of the detector. This effi- 
ciency is estimated to be about 72% near mid-rapidity in 
central events. However, this analysis is found to be in- 
sensitive to variations in the tracking efficiency, as found 
previously [16]. Fake tracks from random combinations 
of hits are generally negligible, reaching only 0.1% for 
|?y| < 1 for the highest multiplicity events. This rate 
increases slightly at large r/. 



IV. DATA ANALYSIS 

Equation 1 implies that each harmonic component of 
the final-state momentum distribution is represented by 
its magnitude v n and azimuthal direction $„. In general, 
any distribution can be expanded into a set of Fourier 
components. However, the distinguishing feature of cor- 
relation due to initial geometry, as opposed to other 
sources of correlations, is that it is a "global" correla- 
tion. That is, $„ specifies a common direction, indepen- 
dent of the particle species, pt and 77, and it drops out 
in the two-particle correlations (Eq. 3). This feature is 
quite different from the correlations expected from jet 
fragmentation or resonance decays, which typically in- 
volve a subset of particles correlated over a finite range 
in A77 with no preferred global direction. Thus, v n can 
be measured either by correlating tracks with the $„ es- 
timated in the forward direction, or it can be measured 
from two-particle correlations with a large A77 gap. In the 
following, the details of these two methods are discussed. 



A. Event plane method 



The azimuthal angle $„ and the beam direction define 
the n th -order reaction plane 2 . However, due to incom- 
plete detector acceptance and finite event multiplicity, 
the true reaction plane angle $ n can not be determined. 
Instead it is approximated by the event plane angle ^„, 
which is defined as the direction of the "flow vector" Q n , 
calculated in this analysis from the Et deposited in the 
FCal towers in each event: 



( ^2 Et cos ncj) — (J2 Et cos n<j>) 
Y^ Et sin rup — (J] Et sin n(j>) 



tann'JJ, 



(5) 



where the sum ranges over towers in the first two lay- 
ers of the FCal (sec Section II). Subtraction of the 
event-averaged centroid removes biases due to detector 
effects [62]. A standard flattening technique is then 
used to remove the residual non-uniformities in the event 
plane angular distribution [63]. These calibration pro- 
cedures are similar to those used by the RHIC experi- 
ments [62, 64]. 

The coefficient v n is measured by correlating tracks 
with \fr„ to obtain the raw values v^° s = (cos n((j> — \fj n )) ■ 



2 If the shape of Pb nuclei is approximated by the smooth Woods- 
Saxon function without fluctuations, <!>„ coincides with the az- 
imuthal angle of the reaction plane defined by the beam axis and 
the impact parameter (the vector separating the barycenters of 
the two nuclei). 



The value of i>° bs is then corrected by a resolution factor 
that accounts for the dispersion of \I/ n about 3> n [4] : 
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where the average is performed over all events for the 
denominator and all tracks and all events for the numer- 
ator. The EP resolution of the FCal, Res{n*„}, is [4]: 



Res{n*„} = (cosn(* n - $„)) 

XnVn -. 



Io(?f) + Ii(^) 



(7) 



where I a are the modified Bcsscl functions of the first 
kind, and Xn (known as the "resolution parameter") is 
the fundamental variable that quantifies the precision of a 
detector for determining the event plane. The value of \n 
is proportional to the -&r-weighted harmonic coefficient 
f^ Cal and the square-root of the total multiplicity M in 
the FCal acceptance [4]: 



Xn OC T^al 



(8) 



The values of Xn and Resjn'I'n} arc obtained from 
a two-subevents method (2SE) and a three-subevents 
method (3SE) [4]. In the 2SE method, the signal of a 
detector used to measure the event plane is divided into 
two "sub-events" covering equal pscudorapidity ranges in 
opposite hemispheres, such that the two subevents nom- 
inally have the same resolution. The FCal detectors 
located at positive and negative 77, FCalp and FCalN, 
provide such a division. The resolution of each FCal 
subevent is calculated directly from the correlation be- 
tween the two subevents: 

Res{n< (N) } = (cosnK (N) - $„)) = v/(coBn(*£ - *$JflB) 

The resolution parameter of the FCal subevent x„ !S ub is 
determined by inverting Eq. 7. The resolution parame- 
ter for the full FCal is Xn = V^Xn,sub, with V2 account- 
ing for a factor of two increase in the total multiplicity 
(Eq. 8). Finally, Xn is incorporated into Eq. 7 to obtain 
the resolution for the full FCal. 

In the 3SE method, the Rcsjn'I'n} value for a given 
subevent A is determined from its correlations with two 
subevents B and C covering different regions in 77: 



Res{nf^} = 



(cosn (tf £ _ tfB)) ( cosn ($A _ ^C)) 



(oosn (*«-*£)) 



(10) 



The large 77 coverage of the ID and ECal, with their fine 
segmentation, allows for many choices for subevents B 
and C. The ID and ECal are divided into a set of 22 
reference subevents each covering 0.5 units in 77. The 
subevents B and C are chosen to ensure a minimum sep- 
aration in 77 of 1 unit between all three subevents. This 
separation in 77 is required to suppress short range cor- 
relations [62]. Various 3SE combinations arc studied to 



check the sensitivity to the size of the chosen pseudora- 
pidity gaps, as well as potential systematic effects due 
to the explicit use of the correlation between FCalp and 
FCal N in the 2SE method. 

Figure 1 shows the values of \ n and Res{n^/„} mea- 
sured as a function of centrality for n = 2-6 using the full 
FCal. The data points and associated statistical uncer- 
tainties are calculated using the 2SE method. However, 
5% upward and 15% downward centrality-independent 
corrections are applied to n = 5 and n = 6 respectively, 
to adjust to the average of the 2SE and the 3SE esti- 
mates for the Rosing}. The differences between the 
two estimates are quoted as systematic uncertainties for 
Rcs{n\I/ n }, and they are propagated via Eq. 7 to obtain a 
systematic uncertainty for Xn- In this analysis, the cen- 
trality range for each harmonic n is chosen such that the 
relative statistical uncertainty for Resjn'I'n} is less than 
30% of its mean value, and the 2SE and 3SE estimations 
show good agreement. They are 0-80% for V2, 0-70% for 
V3 and W4, 0-50% for v§ and vq as indicated in Fig. 1. 

In the event plane analysis, two complementary meth- 
ods are employed to measure v n . The first ("full FCal") 
method calculates w° bs by correlating tracks in the ID 
with the EP from the full FCal detector; the resolution 
correction for the full FCal is then applied to obtain the 
final v n (Eq. 6). In the second ("FCal subevcnt" or sim- 
ply FCalp(N)) method [16], tracks with r\ > (rj < 0) 
are correlated with the EP in the opposite hemisphere 
given by the FCaiN (FCalp). The resolution correction 
for the FCalp (n) is then applied to obtain the final v n . 
Note that the relative statistical and systematic uncer- 
tainties for Res{n 1 I' n } are almost identical for the two 
methods since they both rely on similar subevent cor- 
relations (Eqs. 9 and 10). However, the values of u° bs 
from the FCal subevent method are smaller than those 
from the full FCal method due to its poorer EP resolu- 
tion. Therefore, these v° bs values have a larger fractional 
statistical uncertainty. The primary advantage of the 
FCal subevent method is that it increases the minimum 
(maximum) pseudorapidity separation between the track 
and the EP from about 0.8 (4.8) units for the full FCal 
method to about 3.3 (7.3) units. Thus the subevent ap- 
proach is less affected by short range autocorrelations, 
stemming primarily from jet fragmentation and reso- 
nance decays. In this analysis, the FCal P ( N ) method is 
used for the 7/ dependence of v n to minimize short range 
correlations, while the full FCal method is used for the 
Pt and centrality dependence of v„ to optimize the EP 
resolution (see Section V A). However, the potential influ- 
ence of short range correlations on the full FCal method 
is cross-checked with the FCalp(N) method. Good agree- 
ments are always observed for ^-integrated v n , within the 
systematic uncertainties for the two methods. 

The systematic uncertainty in v° bs is determined by 
varying the track quality cuts, comparing data for dif- 
ferent running periods, varying the full centrality range 
by ± 2% according to the uncertainty in the trigger and 
event selections, as well as by determining the value of 
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FIG. 1. (Color online) The \n (top) and EP resolution factor 
(bottom) vs. centrality (smaller value refers to more cen- 
tral events) for n = 2-6, together with the systematic uncer- 
tainty as shaded bands. The EP is measured by both sides 
of the FCal detector (denote by "full FCal"). Note that the 
Resln^n} value can not be greater than one (see Eq. 7), thus 
its systematic uncertainty shrinks as it approaches one. 



(sinn(4>—^f n )}. The study of track quality cuts accounts 
for influences of background contaminations and tracking 
inefficiency. Finite sine terms can arise from detector bi- 
ases that introduce a correlation between the ID and the 
FCal. Their magnitudes relative to the cosine terms are 
included in the uncertainty for v° bs . All these uncertain- 
ties are generally small, except for n — 6. They are also 
quite similar for the full FCal and FCal subevent meth- 
ods, so the larger of the two is quoted as the main uncer- 
tainty. As a cross-check, V2~vq values are also extracted 
using the EP measured either for the three layers of FCal 
individually or for two r\ regions of FCal (3.3 < | ^-7 1 < 4.0 
and 4.0 < \q\ < 4.8). Although these five FCal subevents 
have up to a factor of four difference in their resolution 
corrections, the measured V2~vq all agree to within 2%- 
10%. 

Tables I-V summarize the systematic uncertainties for 
V2~vq in various centrality intervals. The total uncertain- 
ties are calculated as the quadrature sum of all sources 



in these tables. In most cases, they are specified for mul- 
tiple 5% wide centrality intervals. For example, 0-20% 
in Tables I-V refers to four bins: 0-5%, 5-10%, 10-15%, 
and 15-20%; a "5.0-2.0" notation indicates the values of 
relative systematic uncertainty in percentage at the be- 
ginning and at the end of the 0-20% centrality interval. 
Tables IV and V also quote the uncertainty for the 0-1% 
centrality interval, which generally has the same system- 
atic uncertainty as that for the 0-5% centrality interval, 
but with a larger statistical uncertainty. The system- 
atic uncertainties only include those associated with the 
measurements themselves, and no attempt is made to 
disentangle potential contributions from various sources 
of autocorrelations, as their exact origin and quantitative 
effects on v n are not fully understood [5]. Nevertheless, 
these autocorrelations should be largely suppressed by 
the large average r\ gap between the ID and the detector 
used for determining the EP. 



Centrality 


0-20% 


20-50% 


50-70% 


70-80% 


Resolution [%] 


5.0-2.0 


1.0-2.0 


3.0-4.0 


4.0-6.0 


Track selection [%] 


2.0 


0.5 


0.5 


1.0 


Residual sine term[%] 


0.8 


0.6 


0.5 


0.2 


Running periods[%] 


0.2 


0.2 


0.5 


1.0 


Trigger & event sel.[%] 


1.0 


1.0-0.5 


1.0 


1.5 


Total[%] 


5.6-3.2 


1.4-2.3 


3.4-4.2 


4.6-6.4 



Centrality 


0-20% 


20-50% 


50-70% 


Resolution [%] 


4.0 


4.0 


4.4-16.0 


Track selection [%] 


1.0 


1.0-2.0 


4.0 


Residual sine term[%] 


2.0 


2.0 


3.0-5.0 


Running periods[%] 


1.0 


1.5-2.0 


4.0 


Trigger & event sel.[%] 


0.6 


0.7 


1.0-2.0 


Total[%] 


4.9 


4.9-5.4 


7.9-17.5 



TABLE III. Summary of relative systematic uncertainties in 
percentage for v± for both full FCal and FCal P ( N ). 



Centrality 


0-1% 0-20% 20-40% 40-50% 


Resolution [%] 


10.8 10.2 10.2-10.4 11.2-22.4 


Track selection [%] 


1.0 1.0 1.0 2.0 


Residual sine term[%] 


5.0 


Running pcriods[%] 


2.0 2.0 2.0 4.0 


Trigger & event sel.[%] 


1.0 


Total[%] 


12.1 11.6 11.6-12.1 13.0-23.0 



TABLE IV. Summary of relative systematic uncertainties in 
percentage for v$ for both full FCal and FCal P ( N ). 



TABLE I. Summary of relative systematic uncertainties in 
percentage for vi for both full FCal and FCal P ( N ). See text 
for explanation of the arrangement of the uncertainties. 



Centrality [%] 


0-20% 


20-50% 


50-70% 


Resolution [%] 


3.0 


3.0 


3.0-5.6 


Track selection [%] 


2.0 


0.5 


0.5-2.0 


Residual sine term[%] 


1.0 


1.0 


1.5 


Running periods[%] 


0.5 


0.5-1.5 


2.0 


Trigger & event sel.[%] 


0.4 


0.5-1.0 


1.5-3.5 


Total[%] 


3.8 


3.5-3.9 


4.6-7.4 



Centrality 


0-1% 


0-20% 


20-40% 


40-50% 


Resolution [%] 


58 


34- 


-31 


31 


32-38 


Track selection [%] 


10 


Residual sine term[%] 


10 


Running periods[%] 


10 


Trigger & event sel.[%] 


1 


Total[%] 


61 


38- 


-35 


36 


37-42 



TABLE V. Summary of relative systematic uncertainties in 
percentage for v§ for both full FCal and FCal P ( N ). 



TABLE II. Summary of relative systematic uncertainties in 
percentage for «3 for both full FCal and FCal P ( N j. 



B. Two-particle correlation method 

The two-particle correlation function is generally de- 
fined as the ratio of the same-event pair (foreground) 
distribution to the combinatorial pair (background) dis- 



tribution in two-particle phase space (^aj^bj^aj^b): 

d 4 N 



C(^ a >b,f?a,f?b) 



d 2 N ~ d 2 N 



(11) 



p a dr) a d<ji h dri h 



In practice, the correlation function is usually studied as 
a function of relative azimuthal angle (Ac/)) and relative 
pseudorapidity (A77), by averaging pair distributions over 
the detector acceptance: 



C(A</>,Ar,) 



B(A<b,A v ) ' 



(12) 



where 

S(A0,At7) = 



d(j) a dr] a d(j) h dri h S ah 
d(j) a dr) a d(f) h dri h 6 ah 



d A N 



d(j} a dr) a d(j) h dr] h 

d 2 N d 2 N 

d<p a dr) a d(j) h di] h 



(13) 
(.14) 



The <5 a b is a shorthand notation for 8((f> a — <ph — A<j))8{rj a — 
?7b - Arj). 

For an ideal detector, the combinatorial pair distribu- 
tion is uniform in A<f>, and has a nearly triangular shape 
in Arj due to the weak dependence of the single-particle 
distribution on rj [57]. In reality, both same-event and 
combinatorial pair distributions are modulated by detec- 
tor inefficiencies and non-uniformity. These detector ef- 
fects influence the two distributions in the same way so 
they cancel in the ratio. Therefore, B(A(f>, Arj) is often 
referred to as the pair acceptance function [37]. In this 
analysis, B(A<fi, Arj) is estimated from track pairs from 
two events with similar centrality (matched within 5%) 
and z v tx (matched within 1 mm). The two particles in 
the pair are typically selected with different conditions, 
such as different px ranges, pseudorapidities and charge 
signs. In this analysis, charged particles measured by the 
ID with a pair acceptance extending up to \Arj\ = 5 are 
used. 

Figure 2(a) shows the two-dimensional (2-D) correla- 
tion function for pairs from a representative px range 
of 2-3 GeV and 0-5% centrality interval. It reveals 
the characteristic long-range near-side "ridge" and away- 
side "double-hump" structures that were observed pre- 
viously in heavy ion collisions at the RHIC for a sim- 
ilar px range [41-43]. A narrow short range correla- 
tion peak is also seen at (A<tj, Arj) ~ (0,0), presumably 
due to autocorrelations from jet fragmentation and reso- 
nance decays. From the 2-D correlation function, a one- 
dimensional (1-D) A<j) correlation function can be con- 
structed for a given Arj interval: 



C{Acj>) =ix 



/ S{A(j), Ari)dAri 
J B(A<t>, Ar/)dAr/ 



(15) 



The normalization constant A is determined by scaling 
the number of pairs in 2 < |Ar7| < 5 to be the same be- 
tween the foreground (S) and background (B). This nor- 
malization is then applied to other At? intervals. Each 
1-D correlation function is expanded into a Fourier se- 
ries according to Eq. 2, with coefficients v n . n calculated 
directly via a discrete Fourier transformation (DFT) : 



*N 



= (cosnA(f>) = 



Em=i cos(nA</> m )C(A0 m ) 



-iY 



E„ =1 C(A0 m ) 



.(16) 



where n = 1-15, and N = 200 is the number of A<j) bins. 
A small upward relative correction is applied (~ 0.15% 
for 71 = 6 and increasing to 1% for n = 15) to account for 
the finite A0 bin width. Figure 2(b) shows one such 1-D 
correlation function for 2 < |A?/| < 5, overlaid with the 
corresponding contributions from individual v nn com- 



ponents. The shape of the correlation function is well 
described by the sum of the first six v n , n components. 

According to Eq. 4, if the correlations are dominated 
by those arising from asymmetry of the initial geometry 
such as flow, v nj „ should factorize into the product of 
two single-particle harmonic coefficients. This is found 
to be the case for n > 2 at low px for pairs with a large 
Arj gap, but is not true for n = 1 (see Sections VB and 
VC), similar to what was also found in other measure- 
ments [39, 40]. Thus if the two particles are selected 
from the same px interval ("fixed-px" correlations) as 
in Fig. 2, the single-particle v n for n > 2 can be calcu- 
lated as v n = y/v n!n . When v n _ n < 0, v„ is defined as 

v n = -y/\v n , n \ (or v n = v n , n / y/\v n ,n\ in general). This 
calculation is repeated for all 1-D correlation functions 
in each |A?7| slice. The resulting full |Ar?| dependence of 
v n , n and v n are shown in Figs. 2(c) and 2(d). 

The v n>n and v n values are found to vary rapidly for 
I At; I < 1, presumably reflecting the influence of the short 
range correlation at (A</>, A77) ~ (0,0) (Fig. 2 (a)), but 
they decrease much more slowly for larger |A?j|. This 
slow decrease is expected since the single-particle v n also 
decreases very slowly with 77 (see Fig. 3), and the fac- 
torization relation Eq. 4 is valid for the present px range 
(see Section VB). These behaviors suggest that the auto- 
correlations from near-side jet fragmentation and reso- 
nance decays can be largely eliminated by requiring a 
large A77 gap (e.g. |A?j| > 2). 

Each "fixed-px" correlation function provides a refer- 
ence v n for a chosen px range (denoted by superscript 
"a"). Tracks from this px range are then correlated 
with those from a target px range (denoted by super- 
script "b"), and this "mixed-px" correlation is used to 
calculate u„, ra and to obtain the v„ in the target px via 
Eq. 4. Since factorization is expected to be valid for the 
anisotropics driven by the initial geometry, but is bro- 
ken by the presence of autocorrelations among the jet 
fragmentation products, the level of consistency between 
v n obtained from different reference px ranges reveals 
whether the 2PC is dominated by anisotropics driven by 
the initial geometry. A detailed study of the factorization 
properties of v\-vq is presented in Section VB. 

The correlation function relies on the pair acceptance 
function to reproduce and cancel the detector acceptance 
effects in the foreground distribution. Mathematically, 
the pair acceptance function in A</> is simply a convo- 
lution of two single-particle azimuthal distributions, and 
should be uniform in A<p without detector imperfections. 
A natural way of quantifying the influence of detector ef- 
fects on v ntn and v n is to transform the single-particle 
and pair acceptance functions into the Fourier space. 
The resulting coefficients for pair acceptance v^ c ^ are the 
product of those for the two single-particle acceptances 



,.dct,; 



and v 



dct.b 



In general, the pair acceptance func- 



tion is quite flat: the maximum variation from its aver- 
age is observed to be less than 0.001 for pairs integrated 
over 2 < \Ar)\ < 5, and the corresponding |w^°*| values 
are found to be less than 1.5 x 10 -4 . These v^°* values 
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FIG. 2. (Color online) The steps involved in 
the extraction of v„ values for 2-3 GeV fixed- 
Pt correlations in the 0-5% centrality in- 
terval: (a) two-dimensional correlation func- 
tion, (b) the one-dimensional A<f> correlation 
function for 2 < |Ary| < 5 (re-binned into 100 
bins), overlaid with contributions from the 
individual v n , n components and their sum, 
as well as the residual difference between the 
data and the sum, (c) Fourier coefficient «„,„ 
vs. |Ajj| for n = 1-6, and (d) v n vs. |A»j| 
for n = 2-6. The shaded bands in (c) and 
(d) indicate the systematic uncertainties as 
described in the text. 
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are expected to mostly cancel in the correlation function, 
and only a small fraction contributes to the uncertainties 
of the pair acceptance function. Three possible residual 
effects for u„* are studied: 1) the time dependence of 
the pair acceptance, 2) the effect of imperfect centrality 
matching, and 3) the effect of imperfect z v tx matching. 
In each case, the residual v^% values arc evaluated by 
a Fourier expansion of the ratio of the pair acceptances 
before and after the variation. Overall, significant devi- 
ations are observed only for the effect of imperfect z vtx 
matching, and they are generally larger for narrower |Ar/| 
ranges and higher p^ • 

The systematic uncertainty of the pair acceptance is 
the quadrature sum of these three estimates, which is 
Sv n , n = (2.5-8) x 10~ 6 depending on n, pt, and the 
width of | A?y| interval. This absolute uncertainty is prop- 
agated to the uncertainty in v n , and it is the dominant 
uncertainty when v n is small, e.g. for vq. Moreover, re- 
sults for inclusive charged particles are compared to those 
obtained independently using same-charge and opposite- 
charge pairs. These two types of correlations have some- 
what different pair acceptances due to different relative 
bending directions between the two tracks. They are 
found to give consistent results for n < 6, where the v n ^ n 
values are dominated by physics effects. However, small 
systematic deviations are observed for n > 8, where the 



v n>n values are expected to be dominated by acceptance 
effects. Therefore, the systematic uncertainty also in- 
cludes the RMS difference of the v n , n values averaged for 
8 < n < 15 between the two types of correlations. This 
uncertainty is usually much smaller than those associated 
with v^ni except for large px- 

The second type of systematic uncertainty includes the 
sensitivity of the analysis to track quality cuts, variation 
between different running periods, trigger and event se- 
lection, as well as the ability to reproduce the input v n 
in fully simulated, digitized and reconstructed HIJING 
events with azimuthal anisotropy imposed on the gen- 
erated particles. Most systematic uncertainties cancel 
for the correlation function when dividing the foreground 
and background distributions. The estimated residual ef- 
fects are summarized in Table VI. 

The total systematic uncertainties are the quadra- 
ture sum of the uncertainties calculated from pair ac- 
ceptance, the RMS difference of the i> n , n averaged for 
8 < n < 15 between the same-charge and opposite-charge 
correlations, and those listed in Table VI. They are then 
propagated to uncertainties for v n . These uncertainties 
are plotted as shaded bands around the data points in 
Figs. 2(c) and 2(d). Most of these uncertainties are cor- 
related between different pt ranges. However, a fraction 
of them are found to be uncorrelated with px, coming 





Vl V2 


Vi 


V4 


Vi 


v n\n>6 


Track selection [%] 


3.0 0.3 


0.3 


1.0 


2.0 


4 


Running periods[%] 


5.0 0.3-1.0 


0.7-2.1 


1.2-3.1 


2.3 


7-11 


Trigger & event sel.[%] 


1.0 0.5-1.0 


0.5-1.5 


0.5-1 


1.0 


5 


MC consistency [%] 


2.0 1.0 


1.5 


2.0 


3.5 


5 


Sum of above [%] 


6.3 1.2-1.8 


1.8-3.0 


2.6-3.9 


4.8 


11-14 



TABLE VI. Relative systematic uncertainties for v n in per- 
centage from tracking cuts, variation between different run- 
ning periods, centrality variation, consistency between truth 
and reconstructed v n in HIJING simulation, and the quadra- 
ture sum of individual terms. 



mainly from the track selection, running period varia- 
tion and MC comparison in Table VI and the pair accep- 
tance. This fraction (point to point in pr) is estimated 



to be about 30% of the final systematic uncertainty, and 
the remaining uncertainty is treated as a pr-correlated 
systematic uncertainty. They are used in the discussion 
of the vi i results in Section VD. 



V. RESULTS 
A. v^—ve from the event plane method 

Figure 3 shows the rj dependence of v n for several 
centrality intervals in the 2-3 GeV pt range from the 
FCalp(N) EP method. Similar behaviors are observed in 
other pt ranges (see also [16] for i> 2 ). The t>2 values de- 
crease by less than 5% towards large I77I for central and 
mid-central events, and the decrease is more pronounced 
both for n > 3 and for peripheral events. 

Figure 4 shows the pt dependence of V2~vq for several 
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FIG. 3. (Color online) v n vs. r\ for 2 < px < 3 GeV from 
the FCalp( N ) method (i.e the EP is measured by either FCalN 
or FCalp) with each panel representing one centrality inter- 
val. The shaded bands indicate systematic uncertainties from 
Tables I-V. 



0.2 

0.1 



0.2 
0.1 




0.2 



0.1 



0.2 



0.1 



0-5% ATLAS -*-n=2 

Pb-Pb 1^=2.76 TeV ~°~ n = 3 
-B- n=4 ■ 
Ui= 8 nb" 1 h|<2.5 ^n=5 
full FCal EP -s~n=6 



-+- 



-+- 



10-20% 



30-40% 




H- 



-+- 



50-60% 






5-10% 



9 § 



-+- 



-+- 



■■■?■" 



20-30% 



40-50% 




1 1 1 1 1 



60-70% 



/*N. 




2468 10 2468 10 

P T [GeV] 



FIG. 4. (Color online) v n vs. pt for several centrality inter- 
vals. The shaded bands indicate the systematic uncertainties 
from Tables I-V. 
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FIG. 5. (Color online) v n vs. centrality for six pt ranges from 
the full FCal event plane method. The shaded bands indicate 
systematic uncertainties from Tables I-V. 



centrality intervals. All v n increase with px in the range 
up to 3-4 GeV and then decrease. However, they remain 
positive even at the highest measured px, where occa- 
sional fluctuations to negative values do not exceed the 
statistical precision. This turn-over behavior in px was 
also observed at RHIC for v-i [28, 65], and it is associated 
with the transition from anisotropy driven by the collec- 
tive expansion to anisotropy driven by a path-length de- 
pendent jet energy loss [2, 29]. The overall magnitude of 
v n also decreases with increasing n, except in the most 
central events where v$ is the largest. 

Figure 5 shows the centrality dependence of v n for se- 
veral px ranges. The centrality intervals are presented in 
5% or 10% increments, with an additional interval for the 
1% most central events. Going from central to periph- 
eral events (from right to left along the x-axis), v-i first 
increases, reaching a maximum in the 30-50% centrality 
range, and then decreases. The higher-order coefficients 
v-&-v§ show a similar, but much weaker, centrality depen- 
dence, and this behavior is consistent with an anisotropy 
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FIG. 6. (Color online) V n /v 2 vs. px for several centrality 
intervals. The shaded bands indicate the total systematic 
uncertainties. 



related to the fluctuations in the initial geometry [32]. 
For most of the measured centrality range, vi is much 
larger than the other harmonic coefficients. In central 
events, however, i> 3 and/or uj becomes larger than v 2 for 
some px ranges. At high px (> 4 GeV), t>2 increases to- 
wards more peripheral events, presumably reflecting the 
dominance of autocorrelations from di-jets. 

In an ideal hydrodynamics scenario, v n at low px is a 
power-law function of the radial expansion velocity of the 
fluid, leading to the qualitative expectation that u n (px) 
is a power-law function of px [9, 66]. Previous RHIC re- 
sults have shown that V4/V2 ( or equivalcntly u 4 jv 2 ) 
is almost independent of px [48, 49] 3 . Figure 6 shows 
v n A'2 vs. px for various centrality intervals. These 
ratios vary weakly with px except in the 5% most central 
events, suggesting that such a scaling relation largely ac- 
counts for the px dependence. However, the overall mag- 
nitudes of the ratios seem to vary with centrality and also 
vary slightly with n. 

Figure 7 shows the centrality dependence of v n /v 2 
for 2 < px < 3 GeV. Given that the ratios vary weakly 
with px, the results for other px ranges are similar. The 



3 This V4 was measured relative to the $2 instead of the $4 re- 
action plane, and is known as mixed harmonics [4]. It can be 
regarded as a projection of V4 measured in the $4, onto the $2. 
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FIG. 7. (Color online) v n /v 2 vs. centrality for 2 < pt < 3 
GeV. Lines indicate systematic uncertainty bands, calculated 
by assuming that the uncertainties for different v n are inde- 
pendent. 



ratios are almost independent of centrality in mid-central 
and peripheral events, but then increase sharply toward 
more central events, with a total change of almost a factor 
of two over the 0-20% centrality range. In addition, the 
ratios for n = 4-6 are similar to each other, while they 
are systematically higher than those for n = 3. A similar 
centrality dependence was observed for the v^jv\ ratio 
at the RHIC, and was argued to reflect the centrality 
dependence of fluctuations in the initial geometry [49] . 



B. 



v-2— ve from the two-particle correlation method 



Figure 8 shows the evolution of the 2-D correlation 
function with centrality for particles with 2 < px < 3 
GeV. While central events show structures that are long 
range in Arj (the "ridge" and "double-hump"), the more 
peripheral events show a systematic disappearance of 
these long-range structures and the emergence of clear 
jet-related peaks on the away-side. The magnitude of the 
long-range structures, measured as deviation from unity, 
exhibits a characteristic centrality dependence. Figure 8 
shows that the near-side "ridge" (relative to unity) starts 
at about 0.015 in the 1% most central events, increases 
to 0.12 in 30-50% mid-central events, and then decreases 
and disappears in the most peripheral (80-90%) events. 
Since the harmonics for different n all contribute posi- 
tively to the correlation function at Acf> = 0, this non- 
monotonic centrality dependence simply reflects the fact 
that v n for n > 2 all reach their maxima in the 30- 
50% centrality range for this pt selection, as shown in 
Fig. 5. The away-side long-range structure exhibits a 
similar centrality dependence, but is complicated by the 
contribution from the recoil jet, which starts to dominate 
the away-side shape in the 60-90% centrality range. 

As discussed in Section IV B and shown in Fig. 2, 



each 2-D correlation function is projected onto a set of 
1-D Acf> correlation functions in slices of IA77I, and the 
Fourier coefficients v n ,n and v n are calculated from these 
distributions. Examples of such 1-D correlation func- 
tions are shown in Fig. 9 for pairs with 2 < |At/| < 5 
and 3 < pt < 4 GeV, together with individual contri- 
butions from the first six v n , n components. In a sce- 
nario where the Fourier coefficients are dominated by 
anisotropic flow, the value of the correlation function 
at A</> ~ should be larger than its value at A<fr ~ it 
(see Eq. 3). This indeed is the case up to the 40-50% 
centrality interval, but for centralities greater than 50% 
the trend reverses. This reversing of the asymmetry be- 
tween the near- and away-side amplitudes correlates with 
a continuous decrease of i>i,i, which eventually becomes 
negative at around the 40-50% centrality interval (also 
see top panel of Fig. 11). The correlation function in the 
80-90% centrality interval shows that a broad peak from 
the away-side jet predominantly generate a negative v\ t \ 
and a positive ^2,2- Therefore, they tend to cancel each 
other at the near-side but add up at the away-side. This 
behavior suggests that in peripheral collisions and at low 
Pt, the appearance of a large negative v\ t i is a good in- 
dicator for a significant contribution of autocorrelations 
from jets to ^2,2- 

Figure 10 shows Fourier decomposition of the correla- 
tion functions in the 0-10% centrality interval for several 
Pt ranges. Again, a large pseudorapidity gap of |A?/| > 2 
is required to suppress the near-side jet peak and to ex- 
pose the long-range structures. At low pt, the v n ^ n com- 
ponents are mainly driven by these long-range structures. 
However, for p T ,p T > 6 GeV, they are dominated by the 
pronounced away-side jet correlation centered at A</> = w, 
leading to v n . n values with alternating sign: (—1)". Fig- 
ure 11 contrasts the centrality dependence of the v n _ n for 
a low-pT range (top panel) and a high-pT range (bot- 
tom panel). It shows that at low pt, the sign- flipping of 
!)„.„ between even and odd n happens only for periph- 
eral events; at high pt the sign-flipping happens over the 
entire centrality range. More details on the variation of 
the correlation functions and u„„ values as a function of 
p T and p T are shown in Figs. 29-36 in the Appendix. 

As discussed in Section I, a necessary condition for 
the v ntn to reflect anisotropy associated with the initial 
spatial asymmetries is that it should factorize into the 
product of two single-particle v n . A direct way to verify 
the factorization is to check whether the target v n (p\) 
calculated as: 



, b, Vn,n(PT,PT) VU'n.n (g| , gfQ \Vn,n (Pt, Pt) n ,-s 

Wt) = ~ 7-Ti — = 7. 7ZS, ,^ ( 17 ) 



«u(Pt) 



Vn,n \Pt? P*T ) 



is independent of the reference p T . Figure 12 shows one 
such study of n = 1-6 for 1 < p T < 1.5 GeV using 
four different reference p T ranges. The «i(p T ) values cal- 
culated this way clearly change with the choice of p T , 
indicating a breakdown of the factorization over the mea- 
sured pt range. This breakdown is mainly due to con- 
tributions from global momentum conservation for a sys- 
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tcm with finite multiplicity. More detailed discussions of 
the physics behind this breakdown and ways to separate 
different contributions and extract the v\ associated with 
dipole asymmetry from v\ t \ are presented in Section VD. 
Figure 12 shows that the factorization holds for n = 2- 
6 for | At/ | > 1. Despite a factor of three variation in the 
reference v n value for the px ranges from 0.5-1 GeV to 
3-4 GeV (see Fig. 4), the extracted v n (p T ) for n = 2- 
6 is constant within statistical uncertainties for |Ar/| > 
1. This factorization check is repeated for n = 2-6 for 
various ranges of p T , p T and centrality. For correlations 
with the default choice of |Ar/| > 2 used in this paper, 
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factorization is found to hold at the 5%-10% level for 
Pt < 3-4 GeV in the 70% most central events. Further 
studies of this topic are presented in Section VC. 

Figure 13 shows the extracted Fourier coefficients, v n 
vs. n, for several fixed-px correlations with |Ar/| > 2 
and for the 1% most central events. The p T ' ranges 
are restricted to below 4 GeV, where the factorization 
for V2~vq works reasonably well. Significant signals are 
observed for n < 6, while the signals for n > 6 are con- 
sistent with zero within statistical and systematic uncer- 
tainties. Note that in cases where the uncertainty of u„ jn 
is comparable to its own value, the v n uncertainty be- 
comes asymmetric and highly non-Gaussian. A complete 
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compilation of the Fourier coefficients for other centrality 
intervals is included in Figs 22-27 in the Appendix. Re- 
cent model comparisons indicate that these results can 
provide important constraints on the nature of initial ge- 
ometry fluctuations and the shear viscosity of the created 
matter [34]. 



C. Comparison of v?— ve from the two methods 

Figure 14 compares the centrality dependence of v n 
obtained from the 2PC method and the EP method at 
low px- They agree within 5% for V2-V4, 10% for v$ and 
15% for vq over a broad centrality range, well within the 
quoted systematic uncertainties for the two methods. 

As the contribution of the away-side jet increases with 
Ptj the v n values from the 2PC method are expected 
to deviate from the EP results. Figure 15 compares V2~ 
V5 from the fixed-px 2PC method with those from the 
EP method in the 0-10% centrality interval. The two 
methods agree within 5%-15% for ii 2 and 5%-10% for 
^3 _ ^5 for px < 4 GeV. Significant deviations are observed 
for px > 4-5 GeV, presumably due to contributions from 
the away- side jet (see Fig. 10). 

Figure 16 compares the V2-V5 obtained from the mixed- 
Pt 2PC method with those from the EP method in the 
0-10% centrality interval. The layout of this figure is 
similar to Fig. 15, except that the w n (p T ) are calculated 
using Eq. 17 for four reference p T ranges. Since only one 
particle is chosen from high px, the influence of auto- 
correlations from jets is expected to be less than when 
both particles arc chosen from high px- Indeed, when 



p T is below 3 GeV, the deviations from the EP method 
are less than 5%-10% out to much higher p T , when sta- 
tistical precision allows. This behavior suggests that as 
long as the reference p T is low, the factorization relation 
is valid and hence the v n can be measured to high px via 
the 2PC method. It is also noticed that in central events, 
the agreement for V2 is generally worse than that for the 
higher-order v n : the differences for V2 between the two 
methods reach > 20% for p T > 4 GeV and 3 < p T < 4 
GeV, while the differences for the higher-order v n are 
small and show no visible dependence on p T for p T < 4 
GeV. Comparisons for several other centrality intervals 
arc included in Fig. 28 in the Appendix. Good agreement 
is always observed for V2-V5 in non-central events. 

The agreement between the 2PC and EP methods im- 
plies that the structures of the two-particle correlation 
function at low px and large |A?7| mainly reflect collec- 
tive flow. This is verified explicitly by reconstructing the 
correlation function as: 



C(Acp) 



&^(l + 2^ u cosA0 
2l>r a ^ P ' b cosnA,/>), (18) 



n=2 



where 6 2PC and i>i^ c are the average of the correlation 
function and first harmonic coefficient from the 2PC anal- 
ysis, and the remaining coefficients are calculated from 
the v n measured by the FCalp(N) EP method. Figure 17 
compares the measured correlation functions with those 
reconstructed using Eq. 18, and they are found to agree 
well. In general, the correlation functions arc reproduced 
within 5% of their maximum variations for px < 3 GeV, 
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with slightly larger residual differences for 3 < pt < 4 
GcV. The Vi t i term, not measured by the EP method, 
also plays a significant role, but it is not large enough to 
generate the near- and away-side structures in the corre- 
lation function. 

Figure 17 also shows the separate contributions from 
odd harmonics U3, v§ and even harmonics V2,Vi, vq mea- 
sured by the FCalp(N) EP method. The motivation for 
doing this is that at A0 ~ tt, the even harmonics all give 
positive contributions while the odd harmonics all give 
negative contributions. Thus, the relative magnitudes 
of these two contributions control the away-side shape. 
Figure 17 shows that the "double- hump" structure can 
be explained by the interplay between the odd and even 
harmonics. 



D. ui,i from the two-particle correlation method 

and implications for the v\ associated with the 

dipole asymmetry 

The previous results show that the factorization prop- 
erties of ui t i are quite different from those of the higher- 
order coefficients. The main reason for this is that vi t \ is 
strongly influenced by global momentum conservation, 
while all higher-order coefficients conserve momentum 
due to their multi-fold symmetries. One example of mo- 
mentum conservation effects comes from di-jets in pe- 
ripheral events, which tend to give a negative V\ 1 at 
large |Ar?| (see the cos A0 component in the bottom-right 



panel of Fig. 9). The influence of global momentum con- 
servation on the two-particle correlation has been studied 
extensively [12. 51, 52, 67, 68], and was shown to explic- 
itly break the factorization relation Eq. 4 [51, 52]: 



Vi,\{pT,p\,V^Vb) 



Vl(PT,Va)vi(PT,11b) 



PtPt 

M(p 2 T ) 



(19) 



where M and (p\) are the multiplicity and average 
squared transverse momentum for the whole event, re- 
spectively. The negative correction term is a leading- 
order approximation for momentum conservation. This 
approximation is expected to be valid when the correc- 
tion term is much smaller than one. The global momen- 
tum conservation effect is important in peripheral events 
and for high pt particles, but is diluted in central events 
due to the large multiplicity. 

The first term of Eq. 19 represents a contribution from 
i>i, whose dependence on pseudorapidity can be gen- 
erally separated into a rapidity-odd component and a 
rapidity-even component. The rapidity-odd v\ is thought 
to arise from the "sideward" deflection of the colliding 
ions [5], and changes sign going from negative r\ to pos- 
itive 77 [53, 54]. Since pairs with large |Atj| separation 
typically select two particles from opposite hemispheres, 
the rapidity-odd v\ contribution to V\ t \ is negative at 
large |A?/|. However, the rapidity-odd V\ is found to be 
less than 0.005 for \r\\ < 2 at the LHC [69]. The corre- 
sponding contribution to v\ t \ (< 2.5 x 10 -5 ) is negligible 
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compared to the typical v-^^ values seen in the data (see 
below). The rapidity-even vi signal is believed to arise 
from the dipolc asymmetry of the nuclear overlap due to 
fluctuations in the initial geometry [11, 68]. This spatial 
asymmetry results in a dipole anisotropy of the pressure 
gradient, which drives the v\. If the rapidity-even v\ de- 
pends weakly on 7/ [68] , similar to the higher-order v n , its 
contribution to tji,i should be positive and vary weakly 
with At/. In this case, Eq. 19 can be simplified to: 



vi,i(Pt,Pt) ~ vi{Pt)MPt) ~ 



PtPt 
M{p\) 



(20) 



A previous analysis of published RHIC data, supports 



the existence of a significant rapidity-even v\ at lower 
beam energy [12]. Here the properties of the Vi t i data 
are examined in detail, with the goal of understanding 
the relative contributions of v\ and global momentum 
conservation. 

Figure 18 shows the |Ar/| dependence of the v\ t \ ex- 
tracted from two-particle fixed-px correlations for several 
centrality intervals (see discussions around Fig. 2). The 
D1.1 values for |A?j| < 0.5 are always positive, presum- 
ably due to contributions from the near-side jet or other 
short range correlations. The V\^\ values at large |A?/|, 
however, have a more complex dependence on centrality 
and pt- In peripheral events, the vi t \ values are always 
negative at |Ar/| > 1, and their absolute values increase 
up to I At/ I ~ 1.6 and then decrease for larger |Ar/|. This 
trend is consistent with the gradual falloff of the away- 
side jet contribution at large |Ar/| seen in 2-D correlation 
function, e.g. the last panel of Fig. 8. In central and 
mid-central events, the i>i.i values at |Ar/| > 1 are neg- 
ative for p T ,p T < 2 GeV, but they become positive and 
relatively independent of |Ar/| at higher px, consistent 
with the contribution from a rapidity-even V\ that varies 
weakly with 77. 

Given that the | Ar/| dependence of Fig. 18 is generally 
weak at large |Ar/|, v\.\ is integrated over 2 < |Ar/| < 5 
to obtain one value for each p T and p T combination. The 
v i,i(Pt>Pt) functions are then shown in Fig. 19 for var- 
ious centrality intervals. In peripheral events, the v\^\ 
values are always negative and their absolute values in- 
crease nearly linearly with p T and p\ , presumably dom- 
inated by autocorrelations from the away-side jet. In 
more central events, the absolute value of this negative 
i'i.i component is smaller, reflecting a dilution of the mo- 
mentum conservation term by the large event multiplic- 
ity. Furthermore, Fig. 19 clearly suggests that a positive 
Vi t i component sets in for 2 < px ^S 6 GeV. Its mag- 
nitude increases for more central events and eventually 
drives v\^\ into the positive region. The v\ t i values reach 
a maximum at around 4 GeV, similar to those for the 
higher-order harmonic coefficients (see Fig. 4). 

To validate the two-component model of Eq. 20 as a 
reasonable interpretation of the v\ t i data and to extract 
the vi, a least-squares fit of the «i : i data is performed 
for each centrality interval with the following form: 



E 

a.b 



"1,1 



(p T ,p T ) - [^(ptK'Hpt) - cPtP\])' 



(<?) +(<r 2p 



(21) 



where cr^b* an< ^ cr ab' P P denote the statistical and point 
to point systematic uncertainties for "Ui : i(p T ,p T ), respec- 
tively (see Section IV B). The I'f'^px) function is defined 
via a smooth interpolation of its values at m discrete px 
points, vf'^px.Oli^iJ an d these together with the param- 
eter c constitute a total of m+ 1 fit parameters. An inter- 
polation procedure is used because the functional form of 
vj^ipr) is, a priori, unknown, yet it is expected to vary 
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interval. The error bars and shaded bands represent statistical and systematic uncertainties, respectively. The data points for 
the three highest p^ intervals have coarser binning in p^, hence are connected by dashed lines to guide the eye. The data in 
the bottom three panels are scaled down by a factor of ten to fit within the same vertical range. 



smoothly with p^. 

In the default setup, the Vi t i data used in the fit are re- 
stricted to Pt < 6 GeV and p\ < 10 GeV, giving a total 
of 129 data points for each centrality interval. A cubic- 
spline interpolation procedure is used, and the number 
of interpolation points is chosen to be m = 15 (listed in 
Table VII) . To account for the pT-correlated uncertainty 
(<t/?' c , see Section IV B), the least-squares minimiza- 
tion is repeated after varying all the V\^i data points up 
or down by <r??' corr . Since V\ i(PtjPt) * s symmetric with 
respect to p^, and p\, data points for p^ < p\ are cor- 
related with those for p^ > p\. The effect of this cor- 



relation is evaluated by repeating the fit using only v\ } \ 
data for p^ < p\ or i^i data for p^ > p\. The varia- 
tions from the default fit are included in the systematic 
uncertainties. 

Figure 20 shows the fit to the i>i,i data from Fig. 19 
for the 0-5% and 40-50% centrality intervals. The v\ t \ 
data are plotted as a function of p\ for six intervals of 
p^. The seemingly-complex patterns of the V\ t i data 
are well described by the two-component fit across broad 
ranges of p^ and pJj,, with the dot-dashed lines indicating 
the estimated contributions from the global momentum 
conservation. The typical x 2 /DOF of the fit is between 
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FIG. 20. (Color online) Global fit (via Eq. 21) to the vi,i data for the 0-5% centrality interval (top panel) and the (40- 
50%) centrality interval (bottom panel). The fit is performed simultaneously over all Di,i data points in a given centrality 
interval, which are organized as a function of p\ for various p?^ ranges (indicated at the top of each panel), with shaded bars 
indicating the correlated systematic uncertainties. The fit function and the associated systematic uncertainties are indicated by 
the thick-solid lines and surrounding dashed lines, respectively. The dot-dashed lines intercepting each dataset (with negative 
slope) indicate estimated contributions from the momentum conservation component. The lower part of each panel shows the 
difference between the data and fit (solid points), as well as the systematic uncertainties of the fit (dashed lines). 



one and two depending on the centrality, as shown in 
Table VII. The deviations of the data from the fit, as 
shown in the bottom section of each panel, are less than 
10~ 4 for px < 5 GeV. Above that px and in more pe- 
ripheral events the deviations increase, possibly reflect- 
ing the limitation of the leading-order approximation for 
the momentum conservation term in Eq. 20, or the two- 
component assumption in general. In this paper, the 
study is restricted to the 0-50% most central events, for 
which the first term in Eq. 20 is comparable or larger 
than the second term and the statistical uncertainty of 
the fit is not too large. 

The stability of the least-square minimization proce- 
dure is checked for three types of variations: 1) the inter- 



polation function is varied from cubic-spline interpolation 
to a linear interpolation, 2) the number of interpolation 
points is varied from 9 to 21, and 3) the px range of the fit 
is varied from p T ,p T < 5 GeV to p T ,p T < 10 GeV. The 
systematic uncertainties of the minimization procedure 
are calculated as the RMS sum of these three variations. 
They are small in central events but become substantial 
in peripheral events, and are important sources of uncer- 
tainties at intermediate px- The statistical uncertainty 
of the fit dominates at high px, while the fit uncertainty 
from c^b' 00 " dominates at lowpx- The total absolute un- 
certainty of vf H is dvf* = 0.001-0.004 for p T < 3 GeV 
and increases rapidly for higher px due to the greater 
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FIG. 21. v 1 lt vs. pt for various centrality intervals. The shaded bands indicate the total uncertainty. The uncertainty bands 
are reproduced on their own at the bottom of the figure for clarity. 



statistical uncertainty. 

Figure 21 shows v^ipr) for various centrality inter- 
vals. A significant vf lt; signal is observed for all cases. It 
reaches a maximum between 4 GeV and 5 GeV and then 
falls at higher px ■ This falloff may indicate the onset of 
path-length dependent jet energy loss, which correlates 
with the dipolc asymmetry in the initial geometry similar 
to higher-order v n . The magnitude of the uf 1 * is large: 
its peak value is comparable to that for the V3 shown 
in Fig. 4, and the peak value increases by about 20% 
over the measured centrality range. These results im- 
ply that the rapidity-even collective v\ is an important 
component of the two-particle correlation at intermediate 
Pt- For example, the large positive i>i.i harmonic seen 
in the bottom-right panel of Fig. 17 is mainly due to col- 
lective v\: its contribution to t>i.i is about three times 
larger than the negative momentum conservation term 
estimated by the global fit (top panel of Fig. 20). wf lt; (pT) 
is negative for pt < 1.0 GeV, confirming a generic feature 
expected for collective v\ as suggested by hydrodynamic 
model calculations [11, 68]. This behavior, together with 
the fact that V\\ data show little |A7y| dependence for 
|A7?| > 2 (Fig. 18), is consistent with a rapidity-even v\ 
that is almost independent of rj. 

If the two-component ansatz of Eq. 20 is valid, the fit 
parameter c should be inversely proportional to multiplic- 
ity M and (p T ) of the whole event. This ansatz is checked 
by calculating the product of c and the charged hadron 
multiplicity at mid-rapidity £P- from [57], with the 

assumption that ^p _ is proportional to M. The re- 
sults arc summarized in Table VII for each centrality in- 
terval. Since (p T ) for the whole event is expected to vary 
weakly with centrality (the (px) for charged pions at mid- 



rapidity only varies by ~ 5% within the 0-50% centrality 
interval at the LHC [70]), the product is also expected to 
vary weakly with centrality. Table VII shows that this is 
indeed the case, supporting the assumptions underlying 
Eq. 20. 



Centrality 


X 2 /DOF 


c[x0.001GcV~ 2 


] C ^L [xGcV- 2 ] 


0-5% 


159/113 


0.24 ±0.02 


0.39 ±0.04 


5-10% 


133/113 


0.28 ±0.02 


0.37 ±0.04 


10-20% 


165/113 


0.35 ±0.03 


0.36 ±0.04 


20-30% 


134/113 


0.50 ±0.04 


0.34 ±0.03 


30-40% 


188/113 


0.75 ±0.05 


0.33 ±0.03 


40-50% 


181/113 


1.16 ±0.09 


0.32 ±0.03 


15 interpolation points used in 
0.5, 0.7, 0.9, 1.1, 1.3, 1.5 
3.0, 3.5, 4.5, 5.5, 6.5, 7.5, 


the default fit: 
2.0, 2.5 
9.0 GeV 



TABLE VII. Quality of the fit x 2 /DOF, fit parameter c, and 
corresponding multiplicity scaled values c^- for various 

centrality intervals. The uncertainty of c^P- is calculated 

as the quadrature sum of uncertainties from c and ^P- of 

an 1 17=0 

Ref. [57]. The bottom row lists the 15 pt interpolation points 
used in the default fit. 
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VI. CONCLUSION 



In summary, differential measurements of harmonic co- 
efficients i>i-i>6 for the azimuthal distributions of charged 
particles are presented for lead-lead collisions at ^snn = 
2.76 TeV, based on an integrated luminosity of approx- 
imately 8 /ib -1 recorded by ATLAS. The v n values are 
measured for n = 2-6 using an event plane method, and 
are also derived for n = 1-6 from the v n ^ n measured in a 
two-particle correlation method. 

In the event plane method, V2~vq are extracted as a 
function of transverse momentum (0.5 < px < 20 GeV), 
pseudorapidity (|?7| < 2.5) and centrality (0-80%) by cor- 
relating tracks with the event plane determined at large 
rapidity. The v n values exhibit a weak r\ dependence, 
slightly decreasing toward larger |rj|. The v n values show 
a similar px dependence, namely, increasing with px to 
a maximum around 3-4 GeV and then decreasing for 
higher px. The higher-order coefficients exhibit stronger 
px variations. They follow an approximate scaling rela- 
tion, v n (pt) ot v 2 (Pt)i except in top 5% most central 
collisions. Furthermore, the coefficients v^-vq show little 
dependence on centrality, consistent with an anisotropy 
primarily associated with fluctuations in the initial ge- 
ometry. In contrast, V2 varies strongly with centrality, 
reflecting an anisotropy mainly associated with a change 
in the average elliptic geometry with centrality. 

In the two-particle azimuthal correlation method, har- 
monic coefficients for the distributions of the charged par- 
ticle pairs in A<f>, v n , n = (cosnAcj)), arc extracted over a 
broad range of px, relative pseudorapidity (\Ar]\ < 5) 
and centrality. For pairs of charged particles with a 
large pseudorapidity gap (|Ar/| > 2) and one parti- 
cle at px < 3 GeV, the t>2, 2-^6,6 ar e found to factor- 
ize into the product of corresponding v n harmonics, i.e. 
Vn,n(PT'Pr) ~ v n (p T )v n (p T ), in central and mid-central 
events. This suggests that these values of V2, 2~^6, 6 are 
consistent with the response of the created matter to 
fluctuations in the geometry of the initial state. This 
factorization does not describe the ^2,2^6,6 data in more 
peripheral collisions or at higher px, primarily due to 
autocorrelations from di-jcts. The factorization relation 
is also found to fail for ix.i over the entire px range. A 
detailed investigation of the wi ; i(|A77|,p T ,p T ) data sug- 
gests that they are consistent with the combined contri- 
butions from a rapidity-even v\ component, and a global 
momentum conservation component that increases lin- 
early with the p T x p T and is inversely proportional to 
the event multiplicity. Motivated by this observation, 
a two-component fit is used to extract v\. The derived 
signal is negative at low px, changes sign at px ~ 1.0 
GeV, reaches a maximum at around 4-5 GeV and then 
decreases at higher px- The magnitude of v\ at the peak 
is comparable to that of V3 . 

At low px where single-particle anisotropy is domi- 



nated by collective flow, the two-particle correlation func- 
tions are found to be well reproduced by combining the 
contributions of V2~Vq from the event plane method and 
the wix (Eq. 18). This implies that the main structures 
of the lowpx two-particle correlation for \Ar}\ > 2 largely 
reflect higher-order collective flow (n> 2), together with 
a cosAip term that contains contributions from both a 
rapidity-even v\ and the global momentum conservation. 
Therefore, the low px correlation functions do not al- 
low significant contributions from jet-induced medium 
responses. The fluctuations in the initial geometry to- 
gether with the low viscosity of the created matter are 
potentially responsible for these v n , n coefficients. A de- 
tailed comparison of the comprehensive v n and v n>n data 
presented in this paper with viscous hydrodynamic model 
calculations at low px [25, 26, 32, 34, 71, 72] and jet en- 
ergy loss model calculations at high px [27, 29, 73, 74] 
may help elucidate the nature of these fluctuations and 
better constrain the transport properties of the medium. 
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Appendix A: Comprehensive data plots from 
two-particle correlation 

Section V B only presents the results from selected cen- 
trality or px intervals. One of the primary goals of this 
paper is the detailed study in a wide range of ccntrality 
and pt of the v„,n and v n , obtained from the 2PC in 
A(f> with a large pseudorapidity gap. For completeness, 
the v n coefficients for several fixcd-px correlations and 
for centrality intervals other than that shown in Fig. 13 
are presented in Figs. 22-27. 

One important study of the factorization properties 
of V2~vq is the comparison between the EP method and 
mixed-pT 2PC method. This comparison is shown in 
Fig. 16 for the 0-10% centrality interval. Figure 28 
extends this comparison to several centrality intervals. 
Good consistency is observed for V2~vq in non-central col- 
lisions. 

Figures 29-31 present the full set of the A<fi correlation 



functions with 2 < |Ar/| < 5 for one central, one mid- 
central and one peripheral centrality intervals. These 
correlation functions are the inputs for the v ntn and v n 
spectra. They provide the complete picture of how the 
structures of the correlation functions are influenced, as 
a function of p^ and p\, by the anisotropics related to 
the initial geometry and autocorrelations from the away- 
side jets, as well as how these influences map onto the 
strength and sign of the v n . n coefficients. Note that the 
correlation functions are symmetric with respect to p^ 
and p\ . For example the correlation for 1 < p^ < 2 GeV 
and 2 < p\ < 3 GeV (denoted as [1-2,2-3] GeV) is the 
same as that for 2 < p^ < 3 GeV and 1 < p\ < 2 GeV 
(denoted as [2-3,1-2] GeV). 

Figure 32-36 shows the dependence of t>2,2~V6.6 on p^ 
and p\. Together with Fig. 19, they provide a concise 
summary of v n , n spectra for all the 1-D correlation func- 
tions such as those in Fig. 29-31 for 2 < |A?y| < 5. Since 
the sum of the vi^—vq^ components almost exhaust the 
shape of the correlation functions, these data provide al- 
most equivalent information as the correlation functions 
themselves over broad ranges of pf^ and p\ . 
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